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In this paper we study quantum star graphs with time-dependent bond lengths. Quantum dy- 
namics is treated by solving Schrodinger equation with time-dependent boundary conditions given 
on graphs. Time-dependence of the average kinetic energy is analyzed. Space-time evolution of the 
Gaussian wave packet is treated for harmonically breathing star graph. 

Quantum particle dynamics in nanoscale networks and discrete structures is of fundamental and practical impor- 
tance. Usually such systems are modeled by so-called quantum graphs, systems attracting much attention in physics 
[lj- [3(1 and mathematics during past two decades. 

In physics quantum graphs were introduced as a toy model for studies of quantum chaos by Kottos and Smilansky 
However, the idea for studying of a system confined to a graph dates back to Pauling [4] who suggested to use 
such systems for modeling free electron motion in organic molecules. During last two decades quantum graphs found 
numerous applications in modeling different discrete structures and networks in nanoscale and mesoscopic physics(e.g., 
see reviews (H-jst and references therein). 

Mathematical properties of the Schrodinger operators on graphs inverse problems for quantum graphs [8, 9] 
were also subject for extensive research recently. Also, an experimental realization of quantum graphs was discussed 
earlier in the Ref. [10]. Despite the certain progress made in the study of quantum graphs some of important aspects 
are still remaining as less- or not explored. Especially, this concerns the problems of driven graphs, i.e. graphs 
perturbed by time-dependent external forces. An important example of such a driving force is that caused by driven 
(moving) boundaries. Treatment of such system can be reduced to the Schrodinger equation with time-dependent 
boundary conditions. Earlier, the problem of time-dependent boundary conditions has attracted much attention in 
the context of quantum Fermi acceleration [l^-llj though different aspects of the problem was treated by many 
authors [l6[- [s^l • Detailed study of the problem can be found in series of papers by Makowski and co-authors 21]-(23j. 



It was pointed out in the above Refs. that the problem of ID box with the moving wall can be mapped onto that of 



time-dependent harmonic oscillator confined inside the static box [21]. 

In this paper we treat similar problem for quantum star graph, i.e. we study the problem of quantum graphs 
with time-dependent bonds. In particular, we consider harmonically breathing quantum star graphs, the cases of 
contracting and expanding graphs. The latter can be solved exactly analytically. Motivation for the study of time- 
dependent graphs comes from such practically important problems as quantum Fermi acceleration in nanoscale network 
structures, tunable particle transport in quantum wire networks, molecular wires, different lattices and discrete 
structures. In particular, sites, vertices, nodes of such discrete structures can fluctuate that makes them time- 
dependent. We will study time-dependence of the average kinetic energy and wave packet dynamics in harmonically 
breathing graphs. 

Graphs are the systems consisting of bonds which are connected at the vertices. The bonds are connected according 
to a rule which is called topology of a graph. Topology of a graph is given by in terms of so-called adjacency matrix 

J 1 if i and j are connected . . 

otherwise 

Quantum dynamics of a particle in a graph is described in terms of one-dimensional Schrodinger equation [l|, Q (in 
the units h = 2m = 1): 



dx 2 



k 2 * b (x), b=(i,j), (1) 



2 



where ^denotes a bond connecting zth and jthe vertices, and for each bond 6, the component ^5 of the total wave- 
function ^ is a solution of the eq.(pQ). 

The wavefunction, ^5, satisfies boundary conditions at the vertices, which ensure continuity and current conser- 
vation General scheme for finding of eigenfunctions and eigenvalues for such boundary conditions can be found 
in the Ref.[l|. Different types of boundary conditions for the Schrodinger equation on graphs are discussed in the 
Refs.fsl-Q- I n the following we restrict our consideration by most simplest graph, so-called star graph. The star graph 
consist of three or more bonds connected at the single vertex which is called branching point (see Fig.l). Other ones 
are called edge vertices. The eigenvalue problem for a star graph with N bonds is given by the following Schrodinger 
equation: 



d 2 

— ^(y) = k 2 ^(y), 0<y<lj, j = l,...,7V. 



Here we consider the following boundary conditions [111 ] : 



<Pl\y=0 — <P2\y=0 ~ 

0l|y=Zi = <h\y=h 
N 

E^'lpo = o. 



VN\y=0, 
(j)N\y=l N = 0, 



(2) 



The eigenvalues can be found by solving the following equation [11] 



N 



ctan(k n lj) — 



where corresponding eigenfunctions are given as [ll| 



sin(k n lj) 



sin(k n (lj - y)) 



sr~\ lj+sin(2k n lj) 
sin 2 (k n lj) 



Time-dependent graph implies that lengths of the bonds of a graph are time- varying, i.e., when Lj is a function of 
time. In this case particle dynamics in graph is described by the following time-dependent Schrodinger equation: 



d 2 

i—Vj(x,t) = -—Vjfat), 0<x< Lj(t), j = 1, N : (3) 



with N being the number of bonds. 

In the following we will consider the boundary conditions given by 




^N\x=L N (t) 



These boundary conditions imply that only edge vertices of the graph are moving while center (branching point) is 
fixed. Furthermore, we assume that Lj(t) is given as Lj(t) = ljL(t), where L(t) is a continuous function and lj are 
the positive constants. Then using the coordinate transformation 

x 

y = wr 
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Eq. (j3j) can be rewritten as 

i§i*j(y,t) = -^^oiy^ + i^y^Av^). O < y < l h j = 1, ..,7V. (4) 

It is clear that the Schrodinger operator in the right hand side of Eq. (j4]) is not Hermitian due to the presence of 
second term. Therefore using the transformation 

we can make it Hermitian as 

9 1 <9 2 LL 

We note that the above transformations of the wave function remain the boundary conditions unchanged. 
Time and coordinate variables in Eq. (j5j) can be separated only in case when L(t) obeys the equation 

L 3 L 

- = —C 2 = const, (6) 
In this case using the substitution (pj(y,t) = <j>j{y) • exp (^—ik 2 J j^fp^J •> we get 



For C^O from Eq. (J6j) we have 



^ + (A - C 2 y 2 )4>j = 0, 2/6(0,/,). (7) 



= V^ 2 + /3i + 7 , C 2 = -^(/3 2 -4a 7 ), (8) 
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and 



L(t) = ^/pt + l, c2 = ^ 2 - (9) 

In both cases exact solutions of Eq.(j5j) can be obtained in terms confluent hypergeometric functions. In particular, 
for the case when time-dependence of L(t) is given by Eq.((9j) fundamental solutions of Eq.© can be written as 



and 

= exp M (i - A \-Cy< 
Therefore the general solution of Eq. (j5j) is given as 

<l>j(y)=A j <l> jil +B j <l> ji 2. (10) 
From the boundary conditions given by Eq.(|2|) we have 

Bj = A, A J =A-a J (k), j = 1, 2, 3, TV, 
where B is an arbitrary constant and 

OLjik) = 7, AC h 2 l j = 1,2,. .,7V. 



(a) 



(b) 





FIG. 1: Time-dependence of the average kinetic energy for harmonically oscillating primary star graph. 



Taking to account the relations 



d<t>j,i{y) 



2/=0 



1, 



equation for finding the eigenvalues, k n of Eq.©: 



y=0 



L , from Eq. ([2]) we obtain the following spectral 



f i ^(|-^i-^ 2 ) _cn 



(11) 



Thus the eigenfunctions of Eq.(j5j) can be written as 

fafakn) = ih'ty^id/J + fcd/)], j = 1,2,..., TV. (12) 

Furthermore, we provide the solution of Eq.© for simplest case L(t) =at J r b 1 which correspond to C = in Eq.©. 
In this case of the eigenvalues of Eq.(j5j) can be written as time-dependence of the wall is given as 

A 



sin(y/k n (lj -y)), j = 1,2, ...,iV, 



(13) 



(14) 



where fc n is the nth positive root of the equation 

N 

^2ctan(l?Vk) = 0. 

and L(t) > 0, 5 is the normalization constant. 

Now let us consider harmonically breathing graph, i.e. the case when time-dependence of L(t) is given as 

L(t) = b + a cos uut 

with uj = t^t being oscillation frequency and T is the oscillation period. It is clear that in this case time and coordinate 
variables in Eq. © cannot be separated. Expanding <p(y,i) in Eq. © in terms of static graphs wave functions as 

Vj {y,t) = Y J C n {t)<t>f\y), (15) 
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and inserting this expansion into Eq.(j5j)we have 



C n (t) =^2 M mnC m (t) 

m 

where 

. k 2 m ,LL V 



The quantity we are interested to compute is the average kinetic energy which is defined as 

f)ih ■ (nr. f \ z 

dx. (16) 



i=l n 



1 dipj(x, t) 



dx 



In Fig. 1 time dependence of the average kinetic energy of the harmonically breathing star graph is presented for 
different values of the breathing frequency and amplitude. As it can be seen from these plots, < E(t) > is almost 
periodic for uj = 0.5 and a = 1, while for uj = 10 and a = 1 such a periodicity is completely broken and energy grows 
in time. For uo = 10 and a = 20 the behavior of < E(t) > demonstrates " quasiperiodic behavior". Appearing of 
periodic behavior in < E(t) > can be explained by synchronization of the motion of particle with the frequency. The 
lack of such synchronization causes breaking of the periodicity of the average energy in time. 

Furthermore, we consider wave packet evolution in harmonically breathing star graph by taking the wave function 
at t = (for the first bond) as the following Gaussian wave packet: 

y 1 (x,0) = $(x) = -^e- ! ^ )1 , (17) 



with a being the width of the packet. For other bonds initial wave function is assumed to be zero, i.e. ^(tf, 0) 
^3(^,0) = 0. Then for the initial values of the functions <p^(y,t) in Eq. ([T5]) we have 

/ • \ • L (0)1/(0) 2 

Correspondingly, the expansion coefficients at t = can be written as 



A ^0 



In calculation of the wave packet evolution we will choose initial condition as the wave packet being on the first 
bond only, while for other two bonds the wave function at t = is taken as zero. In Fig. 2 the time evolution of 
the wave packet is plotted for harmonically breathing primary star graph whose bonds oscillate according to the 
law L(t) = 40 + acosut . The oscillation parameters (frequency and amplitude) are chosen as follows: a) lj = 10, 
a = 20; b) uo = 10, a = 1; c) uo = 0.5, a = 1. Fig. 2d presents wave packet evolution in static (time-independent) star 
graph. At t = a Gaussian packet of the width a and velocity is assumed being in the first bond. As it can be 
seen from these plots, for higher frequencies dispersion of the packet and its transition to other bonds occur more 
faster compared to that for smaller smaller frequencies. Again, an important role plays here possible synchronization 
between the bond edge and wave packets motions. Existence or absence of ssuch synchronization defines how the 
collision of the packet with the bond edges will occur and how extensively it gains or loses its energy. Therefore more 
detailed treatment of the wave packet dynamics in harmonically breathing graphs should be based on the analysis of 
the role of synchronization and its cr iter ions. 

In this paper we have treated time-dependent quantum network by considering expanding and harmonically breath- 
ing quantum star graphs. Edge boundaries are considered to be time-dependent, while branching point is assumed 
to be fixed (static). Time-dependence of the average kinetic energy and space-time evolution of the Gaussian wave 
packet are studied by solving the Schrodinger equation with time-dependent boundary conditions. It is shown that for 



(a) 





(c) 



(d) 





FIG. 2: Time evolution of the Gaussian wave packet given by Eq. (|17[) for the parameters: a) uj 
a = 1; c) uj = 0.5, a = 1; d) Wave packet evolution in static star graph. 



10, a = 20; b) u = 10, 



certain frequencies energy is a periodic function of time, while for others it gan be non-monotonically growing function 
of time. Such a feature can be caused by possible synchronization of of the particles motion and the motions of the 
moving edges of graph bonds. Similar feature can be seen also from the analysis of the wave packet evolution. The 
above study can be useful for the treatment of particle transport in different discrete structures, such as molecular 
and quantum wire networks, networks of carbon nanotubes, crystal lattices, and others nanoscale systems that can 
be modeled by quantum graphs. 
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